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EXECUTIVE  SUMMARY 


The  three  major  tasks  we  addressed  in  this  project  were  to:  (a)  identify  the  issues  for 
demonstration  of  a  polymer-based  spin-precession  magnetic  sensor;  (b)  develop  theoretical  and 
computational  techniques  to  study  the  effect  of  the  interface  between  ferromagnetic  metal  and 
polymers  for  FET  applications;  and  (c)  apply  the  existing  light  propagation  codes  to  study 
transmission  of  1.5  pm  wavelength  light  through  GaAs.  We  have  completed  all  three  tasks,  and 
the  major  results  of  our  18-month  effort  are: 

(1)  Detailed  device  modeling  confirms  that  a  room  temperature  detectivity  of  ~  50  (T/a/Hz  is 
possible  with  an  appropriately  fabricated  polymer  magnetic  sensor. 

(2)  Self-consistent  charge  and  spin  transport  calculations  indicate  that  doping  concentrations 
and  field-dependent  mobility  in  the  polymer  and  Schottky  barrier  height  with  the  contacts 
can  be  optimized  for  enhanced  field  effect  transistor  performance  and  magnetic  sensor 
sensitivity. 

(3)  Lao.ySrojMnCfi  (LSMO)  is  demonstrated  to  be  a  ferromagnetic  metal  at  room  temperature. 

(4)  Focused  ion  beam  (FIB)  lithography  has  been  used  to  fabricate  ultra-narrow  trenches 
required  for  high-perfonnance  magnetic  sensors,  but  the  use  of  Ga+  ions  in  the  FIB 
oxidizes  the  ferromagnetic  surface  and  prevents  spin  or  charge  injection  from  half-metal 
Lao.ySrojMnCb  (LSMO). 

(5)  We  demonstrated  the  synthesis  of  very  low  mobility  but  heavily  doped  polymers  and  used 
them  in  magnetic  sensor  fabrication. 

(6)  We  calculated  two-photon  and  free  carrier  absorptioncoefficients  and  changes  to  the 
refractive  index  used  in  light  propagation  codes  and  concluded  that  GalnAs  alloy  is  a  better 
limiting  material  at  1.5  pm  than  GaAs. 

This  report  discusses  our  results  and  conclusions  on  key  requirements. 
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TECHNICAL  DISCUSSION 


1.  SPIN-PRECESSION  MAGNETIC  SENSOR 

This  report  summarizes  the  work  accomplished  from  October  2006  through  March  2008.  The 
main  objective  of  this  work  was  to  identify  the  issues  related  to  realizing  a  polymer-based 
magnetic  sensor  with  ultra-high  detectivity  at  room  temperature.  Identification  of  the  issues 
required  full-fledged  modeling,  fabrication,  and  testing.  In  an  effort  to  leverage  the  findings,  the 
Office  of  Naval  Research  (ONR)  funded  another  iteration  of  fabrication  and  testing.  To  provide 
continuity  and  completeness,  this  report  comprehensively  summarizes  the  results  achieved  under 
both  programs. 

1.1  Device  Concept  and  Operation 

The  proposed  device  is  a  planar  structure  that  uses  half-metal  Lao.7Sr0.3Mn03  (LSMO)  for  the 
contacts  and  an  organic  material  for  electron  transport.  The  contacts  are  magnetically  poled  to 
have  parallel  magnetization.  The  half-metal  nature  of  the  contact  allows  injection  of  only  one 
kind  of  spin  (which  is  parallel  to  the  magnetization)  into  the  organic  material.  The  constituents  of 
the  organics  are  usually  atoms  with  small  atomic  numbers  (such  as  H,  C,  N,  O,  and  S),  and  the 
spin-orbit  coupling  is  extremely  small,  which  allows  the  injected  spins  to  stay  coherent  for  long 
periods.  In  the  absence  of  any  magnetic  field,  the  injected  electrons  will  retain  the  spin  and  can 
find  a  state  in  the  other  contact.  However,  when  a  magnetic  field  is  applied,  the  electron  spin  will 
precess  and  the  spin  orientation  with  the  contact  magnetization  will  change,  resulting  in  an 
increased  resistance.  Since  a  very  small  magnetic  field  is  required  to  alter  the  electron  spin 
orientation,  the  device  is  predicted  to  have  very  high  sensitivity,  even  when  operated  at  room 
temperature. 

1.2  Objectives 

The  objectives  of  this  seedling  effort  are  to 

(1)  Carry  out  more  detailed  modeling,  based  on  self-consistent  spin  transport  through  the 
polymers,  and  predict  the  highest  achievable  sensitivity  of  organic  magnetic  sensors  at 
room  temperature. 

(2)  Identify  the  technical  challenges  related  to  device  fabrication  and  perfonnance,  such  as 
room  temperature  perfonnance  of  half-metal  and  its  contact  resistance;  organic  material 
synthesis,  deposition  and  wetting;  e-beam  etching  for  nanotrenches  in  half  metals;  and 
injection  efficiency  of  earners  from  half-metal  to  organic  material. 

(3)  Identify  the  sources  of  noise,  evaluate  their  effect  on  low-frequency  and  room  temperature 
noise  performance  of  organic  magnetic  sensors,  and  identify  the  technical  challenges  for 
improving  signal-to-noise  perfonnance. 
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1.3  LSMO  Growth  and  Characterization 


LSMO  layers  were  grown  at  Arizona  State  University  by  Professor  Nathan  Newman  and  his 
workers.  The  growth  procedure  used  the  pulsed  laser  deposition  (PLD)  technique,  with  oxygen 
working  pressure  of  350  mTorr,  substrate  temperature  600°C,  laser  power  450  mJ,  frequency 
2  Hz,  and  substrate-target  distance  ~5  cm.  The  duration  of  the  growth  is  2.5  h  to  obtain  300-nm- 
thick  LSMO.  The  grown  film  was  cooled  down  in  a  chamber  with  oxygen  at  a  pressure  of 
~10  Torr  and  a  cooling  speed  of  ~10  K/min.  No  further  treatment  was  done  on  the  as-grown 
LaSrMn03  films.  To  verify  the  effect  on  the  Curie  temperature  (Tc),  two  more  batches  of 
samples  were  grown  in  nearly  identical  conditions,  with  only  the  substrate  temperature  changed 
by  ±  50C. 

We  characterized  the  composition  of  the  LSMO  samples  using  Rutherford  backscattering 
spectroscopy  (RBS)  and  found  it  to  be  very  close  to  the  target  composition  of  33%  strontium. 
The  RBS  channeling  of  as-grown  films  showed  a  minimum  channeling  parameter  of  ~4%, 
indicating  excellent  epitaxial 
growth  on  the  SrTi03  substrate. 

One  small  sample  in  each  batch 
was  characterized  by  the  physical 
property  measurement  system 
(PPMS)  to  get  the  Curie 
temperature.  The  Tc  was  found  to 
be  nearly  independent  of  the 
substrate  temperature.  The 
measured  magnetic  moment  as  a 
function  of  temperature  is  shown  in 
Figure  1.  We  see  that  Tc  is  close  to 
360  K.  The  theoretical  maximum 
for  high-quality  LSMO  is  about 
370  K. 

1.4  Fabrication  of  Nanotrenches 

Our  preliminary  design  required  trenches  between  50  and  200  mn  thick  in  the  LSMO  layer.  The 
LSMO  on  either  side  of  the  trench  forms  the  half-metal  contact.  Initially  we  explored  e-beam 
lithography  with  appropriate  masking  as  a  means  to  make  narrow  trenches.  Although  e-beam 
lithography  can  in  principle  be  used  for  this  purpose,  it  involves  multiple  and  intricate  processing 
steps.  We  found  a  recently  developed  technique,  focused  ion  beam  (FIB)  lithography,  to  be  more 
cost  and  time  efficient  in  making  narrow  trenches.  We  have  successfully  used  FIB  lithography  to 
make  trenches  varying  in  thickness  from  ~1  pm  (Fig.  2a)  to  ~25  mn  (Fig.  2b).  In  addition,  we 
observed  the  trench  walls  fabricated  using  FIB  to  be  smooth  and  clean,  a  key  requirement  for 
enhanced  sensor  sensitivity. 


Figure  1:  Magnetic  properties  characterization  of  one  small 

sample  in  Batch  1  by  PPMS.  The  Curie  temperature  is 
about  360  K.  Thickness  of  this  film  is  around  300  nm. 
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We  made  a  mask  design  with  40  devices  per  1  cm  x  1  cm  strontium  titanate  (STO)  substrate. 
Each  device  is  approximately  200  pm  wide.  To  electrically  separate  the  devices,  we  etched 
approximately  50-pm-thick  LSMO  between  the  devices.  The  FIB  technique  is  not  appropriate 
for  efficient  removal  of  a  large  volume  of  material.  We  found  that  Ar-ion  plasma  etching  did  not 
work  well  and  the  wet  etching  with  HC1+KI  removed  the  LSMO  very  quickly  and 
uncontrollably.  We  finally  used  accelerated  Ar-ion  etch  and  then  wire  bonded  the  LSMO 
contacts.  Fully  fabricated  devices  were  supplied  to  Johns  Hopkins  University  for  polymer 
deposition  and  magnetic  testing.  In  the  second  round  of  experiments,  SRI  deposited  similar 
polymers  and  conducted  electrical  testing  of  the  devices. 


(a)  (b) 


Figure  2:  Microtrench  (a)  and  nanotrench  (b)  made  on  LSMO  with  FIB  lithography. 

1.5  Polymer  Synthesis 

The  polymer  synthesis,  deposition,  and  characterization  described  in  this  section  were  carried  out 
at  JHU  by  Professor  Howard  Katz  and  coworkers.  Calculations  carried  out  at  SRI  under  the 
DARPA  program  indicated  that  polymers  with  low  mobility  (~10’  cm'/ Vs)  and  high  carrier 
density  (~10  cm'  )  are  optimal  for  enhanced  magnetic  sensitivity  in  the  proposed  spin 
precession  device.  In  a  recent  literature  report  [Majumdar  et  ah,  Appl.  Phys.  Lett.  89,  122114, 
2006],  the  authors  demonstrated  spin  injection  from  LSMO  into  regioregular  poly 
3-hexylthiophene  (P3HT).  The  mobility  and  carrier  density  of  this  polymer  are  about  10'4  cm2/Vs 
and  10  cm'  ,  respectively.  Hence  we  considered  this  polymer  as  a  starting  point  for  modifying 
the  characteristics  as  required  for  our  sensor  studies. 

The  blended  polymers  are  deposited  on  a  Si  FET  structure  to  measure  the  drain  current  as  a 
function  of  drain  and  gate  voltages.  From  the  well-known  formula  for  field-effect  transistor 
currents  in  saturation, 

W=^C„,i„.(Va- VT I2 

and  the  known  values  of  width,  length,  and  capacitance,  we  obtain  the  mobility  of  the  polymer. 


4 


Then,  from  measured  resistance  R  and  from  the  following  formula  we  can  obtain  the  free  carrier 
density. 

*— L' 

nq/j  A 

n  =  free  carrier  density  in  cm  3 

g  =  1.6xl0  19  C 

ju  =  mobility  in  cm2 1  Vs 

d  =  film  thickness  between  source  and  drain  in  cm 
A  =  cross  -  sectional  area  of  the  current  path  in  cm2 
R=  Resistance  in  ohm 

We  carried  out  two  parallel  studies.  In  one  study,  we  systematically  blended  the  insulator 
polypropylene  (PP)  into  P3HT.  Table  1  lists  electrical  properties  obtained  from  films  on  Si/SiCh 
with  gold  electrodes  on  top  and  an  aspect  ratio  of  the  region  between  electrodes  of  about  20. 


Table  1 .  Electrical  properties  of  P3HT-polypropylene  blends  in  transistor  geometries. 
Resistance  is  the  two-terminal  resistance  with  zero  gate  voltage. 


Mobility  (cm2/Vs) 

Resistance  (MQ) 

Charge  Density  (/cm3) 

P3HT  only 

3.2  x  1 0"4 

118 

1.4  x  1017 

PP  1  % 

5.3  x  10'4 

31 

3.1  x  1017 

PP  2% 

6.0  x  10"4 

32 

2.7  x  1017 

PP  5% 

2.1  x  10'4 

160 

1.5  x  1017 

PP  1 0% 

5.1  x  10'4 

66 

1.5  x  1017 

PP  25% 

0.8  x  10'4 

187 

3.4  x  1017 

For  low  concentrations  of  PP  considered  here,  neither  mobility  nor  carrier  concentration  changed 
systematically.  However,  as  the  PP  content  is  increased  to  25%,  the  mobility  decreased  4-fold 
and  carrier  density  increased  2-fold.  With  more  PP,  the  mobility  is  expected  to  decrease  even 
further  without  affecting  the  carrier  density.  Space-charge-limited  current  measurements  in 
vertical  geometries  are  consistent  with  these  conclusions. 

In  the  other  study,  we  systematically  added  HBF4,  an  acid  dopant  intended  to  increase  the  carrier 
density.  Table  2  lists  the  results. 
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Table  2.  Electrical  properties  of  P3HT-HBF4  blends  in  transistor  geometries. 


Mobility  (cm2/Vs) 

Resistance  (MQ) 

Charge  Density  (/cm3) 

HBF4  1% 

2.4  x  10"2 

4 

4.9  x  1016 

HBF4  2% 

0.6  x  10"2 

65 

1 .2  x  1016 

HBF4 5% 

1.7  x  10  '2 

9 

3.3  x  1016 

HBF4  10% 

1.9  x  10  '2 

4 

6.0  x  1016 

The  measured  mobilities  are  nearly  constant  and  higher  than  the  desired  values,  and  the 
measured  charge  densities  show  no  correlation  with  the  dopant  concentration.  However,  by 
adding  PP  to  the  dopant  blends,  the  mobilities  are  expected  to  decrease  and  lower-mobility 
charge  carriers  will  count  toward  the  charge  density.  Thus,  we  expect  that  an  appropriate 
combination  of  PP  and  HBF4  in  P3HT  will  yield  much  lower  mobility  and  much  higher  carrier 
density  than  those  of  bare  P3HT.  This  will  be  considered  in  future  studies. 

We  obtained  even  lower  mobilities  (than  P3HT-HBF4)  with  higher  charge  densities  by 
substituting  an  organic  electron-acceptor  dopant,  tetrafluorotetracyanoquinodimethane 
(F4TCNQ).  Results  are  shown  in  Table  3.  A  combination  with  appropriate  PP  content  will  be 
tested  in  future  synthesis. 


Table  3.  Electrical  properties  of  P3HT-F4TCNQ  blends  in  transistor  geometries. 


Mobility  (cm2/Vs) 

Resistance  (MQ) 

Charge  Density  (/cm3) 

P3HT  only 

4.6  x  10-4 

18 

6.2  x  1017 

F4-TCNQ  2% 

2.7  x  10-3 

5 

3.4  x  1017 

F4-TCNQ  5% 

3.5  x  10-3 

5 

2.6  x  1017 

F4-TCNQ  10% 

4.3  x  10-3 

4 

3.1  x  1017 

We  attempted  to  characterize  the  resistance  of  P3HT  cast  between  LSMO  electrodes,  both  in 
nanogaps  and  between  macroscopic  electrodes.  In  the  absence  of  polymers,  the  device  resistance 
is  in  gigaohms.  With  P3HT  deposited,  there  should  have  been  increase  in  current  and  decrease  in 
resistance.  In  one  case,  deposition  between  two  macroscopic  electrodes  led  to  device  resistance 
of  ~50  MQ,  convincingly  demonstrating  a  charge  injection  through  polymer.  However,  when  the 
polymer  was  deposited  on  devices  with  nanotrenches,  we  did  not  see  any  noticeable  increase  in 
current. 

First  we  speculated  that  PP,  which  is  an  insulator,  got  attached  to  the  LSMO  walls,  thus 
preventing  any  charge  injection.  The  device  was  cleaned  with  chloroform,  which  removed  all 
polymers.  Then  we  drop-cast  P3HT  without  PP  or  HBF4  on  the  devices.  We  did  not  see  any 
change  in  current.  We  then  further  modified  our  deposition  process  to  1)  change  the  solvent  from 
chlorobenzene  to  chlorofonn,  which  has  lower  boiling  point;  2)  use  a  diluted  polymer  solution 
(2  mg/mL)  for  drop-casting;  3)  use  two  kinds  of  RR-P3HT  polymer,  one  with  a  molecular  weight 
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of  87,000  from  Aldrich  and  the  other  with  a  molecular  weight  of  70,000  from  Rieke  Metals,  Inc.; 
4)  deposit  polymer  under  house  vacuum  before  the  evaporation  of  chloroform  mother  solution; 
and  finally  5)  remove  the  residual  solvent  by  annealing  at  85°C  for  over  an  hour.  In  spite  of  these 
changes  to  the  deposition  methods,  we  did  not  see  any  increased  current  appropriate  for  the 
resistance  of  the  polymer  in  the  nanotrench.  We  speculated  that  the  polymers,  deposited  by  spin- 
coat  and  drop-cast,  did  not  get  into  the  nanotrenches,  or  that  the  LSMO  electrodes  retained  an 
insulating  film  from  their  fabrication. 

SRI  followed  up  the  conclusions  from  JHU  with  a  local  effort  to  deposit  commercially  available 
polymer  (P3HT)  with  a  molecular  weight  of  17,500  and  followed  the  deposition  procedure  as 
described  above.  First,  we  deposited  the  polymer  in  the  microgap  between  the  LSMO  lines.  The 
side  walls  were  coated  with  oxide,  and  the  current  from  LSMO  was  expected  only  from  the  top 
exposed  surface.  The  current  voltage  was  measured.  The  resistance  calculated  for  this  geometry 
with  the  measured  values  of  mobility  and  carrier  concentration  matched  up  very  well  (mostly 
within  20%)  of  the  measured  resistance.  This  study  confirmed  the  charge  injection  through 
polymer  from  LSMO.  Second,  we  deposited  the  polymer  in  the  nanotrench.  Our  electrical 
measurements  clearly  indicated  the  presence  of  polymer  in  the  nanotrench.  However,  the 
measured  current-voltage  curve  was  highly  non-ohmic  in  that  the  current  did  not  increase 
linearly  with  the  cross-sectional  area  or  decrease  linearly  with  trench  width.  This  led  us  to 
conclude  that  the  charge  injection  from  LSMO  to  the  polymer  is  through  an  insulating  oxide 
layer(s).  Further  literature  search  revealed  that  LSMO  reacts  with  Ga-ion  to  form  an  insulating 
and  non-ferromagnetic  layer  [Ono  et  al.,  Appl.  Phys.  Lett.  84,  2370,  2004],  Since  the  nanotrench 
was  fabricated  using  the  FIB  method,  which  uses  a  high-energy  Ga-ion  beam,  we  infer  that  an 
insulating  layer  created  at  the  surface  prevented  the  charge  injection.  We  concluded  that  e-beam 
lithography  is  required  to  make  the  nano  trenches  in  LSMO.  In  this  approach,  an  e-beam  is  used 
to  make  an  ultra -narrow  trench  on  an  e-beam  photoresist  deposited  on  LSMO.  Then,  energetic 
AL  ion  beam  milling  (IBM)  is  used  to  etch  the  LSMO  away.  The  important  requirement  for  the 
e-beam  resist  is  that  its  IBM  etch  rate  be  slower  than  that  of  LSMO. 

1.6  Noise  Modeling 

The  noise  measurement  and  modeling  described  in  this  section  were  carried  out  by  Professor 
Nowak  at  the  University  of  Delaware.  In  the  absence  of  devices  showing  both  charge  and  spin 
injection  into  the  polymer,  we  focused  on  characterizing  the  noise  properties  of  LSMO  thin  films 
and  regioregular  P3HT  separately.  These  studies  reveal  some  of  the  noise  mechanisms  relevant 
to  devices  and  hence  provide  a  rudimentary  model  for  the  device  noise.  Additional  noise  sources 
in  actual  devices  may  include  surfaces  of  the  LSMO  sidewalls  (due  to  lithographic  processing) 
and  the  presence  of  an  interfacial  layer  between  the  LSMO  sidewall  and  polymer.  We  note  that 
noise  studies  of  the  polymer  were  done  on  macroscopic  thin  films.  The  noise  of  a  polymer  in  the 
trenches  may  be  different  if  the  width  of  the  trench  is  smaller  than  the  length  scale  characterizing 
the  disorder  (e.g.,  grain  boundaries)  in  the  polymer  film.  Next,  we  describe  the  noise  in  patterned 
LSMO  lines  and  in  P3HT  that  was  spin-cast  between  the  gold-coated  LSMO  contact  pads. 
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SRI  patterned  LSMO  films  lithographic-ally  into  long,  narrow  lines  and  post-deposited  gold 
contact  pads  onto  the  ends  of  the  LSMO  lines.  Electrical  leads  were  ultrasonically  bonded  to  the 
contact  pads.  The  contact  resistance  of  the  bonds  was  ohmic  and  measured  to  be  less  than  3.5  f2 
per  bond  pad,  far  smaller  than  the  LSMO  line  resistance  of  ~18  kO.  Five  LSMO  lines  (6,  7,  8, 
12,  and  an  unknown  line)  were  measured.  Both  2-probe  and  4-probe  techniques  were  used  to 
measure  the  resistance,  current-voltage  characteristics,  and  noise.  The  results  were 
indistinguishable,  indicating  the  contacts  contributed  negligibly  to  the  measured  noise  power. 
Figure  3  shows  a  typical  linear  I-V  characteristic.  The  resistance  of  the  lines  ranges  from 
approximately  17  to  19  kO,  which  corresponds  to  an  LSMO  resistivity  of  7  mfl/cm. 


Figure  3:  Current-voltage  curve  for  line  12. 


Figure  4  shows  the  voltage  noise  power  spectral  density  as  a  function  of  frequency  at  different 
current  biases.  In  this  plot  the  thermal  (Johnson)  noise  of  the  sample  has  been  subtracted.  The 
measured  thermal  noise  power  spectral  density  was  in  excellent  agreement  with  the  prediction, 
4 kBTR .  The  spectral  dependence  shown  in  Fig.  4  is  very  close  to  \/f  at  low  frequencies  (below  1 
kHz).  Also,  the  noise  power  is  found  to  scale  as  the  square  of  the  bias  current,  /,  indicating  the 
noise  is  due  to  resistance  fluctuations  in  the  LSMO  material.  A  Hooge  expression  can  be  defined 
to  parameterize  the  noise: 


V 

Sv  =  a 

Qf 
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Figure  4:  Frequency  dependence  of  the  voltage  noise  power  in  LSMO  line  12 
measured  for  various  bias  currents.  Spikes  near  60  and  120  Hz  are 
artifacts  associated  with  60  Hz  power  supply  lines  coupling  into  the 
measurement  circuit. 


In  this  expression,  Q  is  the  volume  of  the  sample,  Vdc  =  IR  ;  f  is  frequency;  and  the  quantity  a 
parameterizes  the  magnitude  of  the  noise.  The  quantity  a  allows  comparison  of  noise  on  different 
samples  potentially  measured  under  different  bias  conditions.  We  note  that  the  inverse  scaling 
with  sample  volume  is  plausible  but  has  not  been  explicitly  verified  in  our  samples  since  they  all 
had  the  same  volume.  Of  the  five  LSMO  lines  measured,  for  three  of  them 
a  =  2x\0~10 jum3  ±lxlO“10///?z3 ,  whereas  for  line  6  a  =  lx\0~9^im3,  and  for  the  unknown  line 
gave  a  =  1  xl(T8/wr3 .  We  believe  a  =  2x\0~w jum3  is  a  representative  value  for  LSMO.  This  value 
is  also  consistent  with  some  reports  in  the  literature.  Figure  5  shows  the  magnetoresistance  (MR) 
of  LSMO  line  12  for  the  magnetic  field  applied  along  the  length  of  the  line  and  perpendicular  to 
the  length.  The  MR  is  quite  small  in  both  cases,  being  approx-imately  0.4%  for  H  perpendicular 
to  the  line  and  0.8%  for  H  parallel  to  the  line. 


Field  (G) 
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Figure  5:  Magnetoresistance  of  LSMO  line  12  for  magnetic  field  applied 

parallel  and  perpendicular  to  the  length  of  the  line.  Curves  offset 
for  clarity. 

The  magnetic  field  dependence  of  the  noise  power  (measured  in  an  octave  centered  at  1 1  Hz)  is 
shown  in  Figure  6.  The  noise  exhibits  negligible  dependence  on  magnetic  field.  In  principle, 
thennal  magnetization  fluctuations  can  couple  to  the  resistance  (via  [dR/dM )  ,  where  Mis  the 
magnetization)  and  generate  resistance  fluctuations.  Under  quasi-equilibrium  conditions,  the 
signature  feature  of  this  magnetic  noise  is  that  the  resistance  noise  power  varies  as  dR/dH,  the 
derivative  of  the  magneto-resistance  curve  in  Figure  5.  We  have  not  observed  this  characteristic 
feature.  The  most  likely  reason  that  1  If  magnetic  noise  is  not  evident  in  these  LSMO  lines  is  that 
the  magnetoresistance  is  small.  Hence,  we  attribute  the  1  If  noise  in  the  LSMO  as  being  due  to 
scattering  from  defects  which  undergo  thermally  activated  reconfigurations  -  this  is  the 
commonly  held  explanation  for  1  If  noise  in  metals. 
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Figure  6:  Magnetic  field  dependence  of  the  noise  power  measured  near  1 0  Hz  for 
LSMO  line  12.  Data  was  taken  with  a  constant  current  bias  of  125  pA  at 
room  temperature. 

It  is  important  to  note  that  the  following  data  on  P3HT  were  taken  on  a  different  die  than  the  one 
on  which  noise  in  LSMO  was  measured.  The  noise  in  the  P3HT  polymer  was  measured  in 
samples  that  had  the  polymer  either  spin-cast  or  drop-cast  across  several  LSMO  lines.  In 
principle,  the  polymer  was  confined  to  a  small  (but  macroscopic)  window  that  was  opened  up  in 
the  SiOx  or  SiNx  overlayer  that  covers  the  LSMO  lines.  After  examining  the  polymer  under  an 
optical  microscope,  however,  we  were  concerned  that  additional  polymer,  or  some  unknown 
residue,  had  also  bridged  a  number  of  the  contact  pads.  If  so,  the  current  distribution  is  poorly 
defined  and  it  is  difficult  to  draw  conclusions  related  to  where  the  noise  source(s)  are  and  to 
establish  benchmark  parameters  for  the  polymer.  Moreover,  since  LSMO  material  is  used  as  a 
contact  to  the  polymer,  not  all  the  noise  can  be  simply  attributed  to  the  P3HT. 
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Typical  I-V  characteristics  across  different  pairs  of  LSMO  having  P3HT  between  them  are 
shown  in  Figure  7.  The  curves  are  nonlinear.  The  resistance  between  adjacent  LSMO  lines  (lines 
23,  24,  25,  26)  is  of  order  17  MO  and  scales  approximately  linearly.  That  is,  7?=16.7  MO 
between  lines  23-24;  R  =  32.5MO  between  23-25;  and  R  =  47  MO  between  23-26.  The  data 
suggest  that  there  is  charge  injection  into  P3HT. 


The  noise  data  we  took,  however,  were  not  reproducible  and  showed  peculiar,  non-monotonic 
dependence  on  the  bias  current  (see  Figure  8).  Noise  measurements  between  some  adjacent 
LSMO  lines  showed  day-to-day  variations  in  noise  power  as  large  as  a  factor  of  100.  This 
peculiar  behavior  is  not  understood.  It  is  likely  that  the  nonunifonn  current  distribution  in  the 
polymer,  microstructural  aging  an  unintentional  doping  from  the  atmosphere,  have  a  drastic 
effect  on  the  measured  transport  properties.  Consequently,  we  draw  no  conclusions  regarding  the 
noise  in  P3HT  at  the  present  time. 
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Figure  8:  Voltage  noise  power  spectral  density  in  P3HT  measured  across  adjacent 
LSMO  lines  23  and  24. 

The  final  device  will  consist  of  a  series  of  two  LSMO  lines  sandwiching  a  polymer.  The  voltage 
noise  power  in  the  device  can  be  modeled  as  the  superposition  of  independent  noise  sources.  At 
the  moment  we  anticipate  that  the  dominant  sources  of  noise  will  be  thermal  Johnson  noise 
related  to  the  equivalent  resistance  of  the  entire  device,  1  If  resistance  noise  in  the  LSMO,  and  I  If 
noise  in  the  polymer.  We  can  safely  ignore  amplifier  noise  because  the  equivalent  resistance  of 
the  device  will  be  of  order  10  MO,  which  yields  a  thermal  Johnson  noise  level  of  400  nV/vHz, 
well  above  the  noise  level  at  1  Hz  of  a  decent  differential  amplifier.  Putting  it  all  together,  the 
voltage  noise  power  spectral  density  of  the  device  is  given  by: 


'i  device 


Johnson  +  g\!  f,LSMO  +  gllf,P3HT 


s  device  =  +  ^  '  DC , LSMO  +  ^ !  f  ,P3HT 


v 1 

v  n 


Of 

The  equivalent  magnetic  field  noise  power  in  teslas  per  hertz  can  be  computed  using  the 
response  of  the  device  to  magnetic  field  ( dV/dB );  namely: 
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Most  of  the  voltage  drop  across  the  device  will  take  place  across  the  P3HT,  and  so  the  noise 
contribution  from  LSMO  can  be  comparatively  small.  Future  work  will  be  directed  at 
quantifying  the  noise  in  P3HT  in  samples  having  well-defined  geometries  with  uniform  current 
distributions  and  using  low-noise  metallic  contact  pads. 
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1.7  Improved  Spin  and  Charge  Transport  Modeling 

We  had  previously  considered  the  electric  field  inside  the  organics  to  be  constant  and  the  charge 
distribution  to  be  homogeneous.  In  this  approximation,  the  voltage  dropped  uniformly  across  the 
device.  This  resulted  in  separate  linear  differential  equations  for  spin  and  charge  transport.  We 
have  improved  the  calculations  to  remove  those  approximations.  When  the  approximations  are 
removed,  we  find  that  the  required  differential  equations  are  nonlinear,  as  shown  below. 

V.E  =  _— (p„-ND)  (1) 
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Also  notice  that  the  charge  transport  (Eq.  2)  and  the  spin  transport  (Eq.  3)  have  now  to  be  solved 
self-consistently  as  they  are  coupled  through  the  Poisson  equation  (Eq.  1).  We  have  completed 
the  calculations  and  found  that  linear  approximation  is  valid  when  the  organic  has  a  doping 
density  of  6  x  10  cm'  or  more.  For  lower  dopant  densities,  the  field  is  far  from  uniform. 

We  calculated  the  response  function,  dV/dB,  the  quantity  that  directly  relates  to  sensitivity.  The 
results  are  shown  in  Figure  9  for  two  dopant  concentrations:  6  x  10  cm'  (Fig.  9a)  and 
1x10  cm'  (Fig.  9b).  We  see  that  the  previously  estimated  sensitivity  (red  line)  agrees  mostly 
with  the  accurately  evaluated  value  (blue  line)  only  at  high  doping  density.  This  clearly  indicates 
that  highly  doped  organic  is  needed  for  enhanced  sensitivity. 

These  calculations  and  results  have  been  summarized  in  the  form  of  a  manuscript  for  possible 
publication  in  Physical  Review  B.  The  draft  manuscript  is  enclosed  here  as  Appendix  A. 


(a)  (b) 

Figure  9:  Response  function  calculated  with  linear  (blue)  and  nonlinear  (red)  approximation  for  dopant 
density  of  (a)  1  x  1 018  cm'3  and  (b)  1  x  1 018  cm'3. 
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1.8  Performance  Modeling 


The  sensitivity  of  devices  is  limited  by  several  noise  sources  related  to  material  properties, 
device  design,  interface  quality,  volume,  temperature,  concept  of  operation,  instrument,  etc.  The 
noise  sources  fall  into  four  general  categories:  shot  noise,  Johnson  noise,  1/f  noise  with  magnetic 
origin,  and  1/f  noise  with  nonmagnetic  origin.  The  performance  prediction  depends  critically  on 
the  model  used  to  describe  each  noise.  With  additional  funding  to  this  seedling  effort  we  will 
develop  a  database  from  our  measurements  for  noise  modeling.  Until  then  we  will  use  the  noise 
model  developed  for  organics  from  field-effect  transistor  (FET)  measurements  [Martin  et  ah,  J. 
Appl.  Phys.  87,  3381  (2000)]  and  the  spin  diffusion  length  estimated  from  spin  injection 
experiments  [Dediu  et  al,  Solid  State  Comm.  122,  181  (2002)]. 

The  organics  used  in  these  experiments  differ  from  the  one  we  intend  to  use  in  our  magnetic 
sensor  in  three  important  ways.  First,  the  FET  uses  crystalline  organics  with  high  mobility 
(-1  cnr/Vs).  In  our  work,  we  will  use  low-mobility  (—10'  cnr/Vs)  organics  in  amorphous  state. 
The  dominant  scattering  mechanism  in  crystalline  materials  need  not  be  the  limiting  one  in  the 
disordered  state. 


Second,  the  molecule  T(,  used  in  spin  injection  experiments  contains  S  atoms.  A  spin  coherence 
length,  Ls,  of  50  mn  explains  their  results.  We  intend  to  use  polymer  MEH-PPV  in  which  the 
heaviest  atom  is  O.  Because  of  the  smaller  atomic  number,  the  spin-orbit  coupling  is  expected  to 
be  weaker  and,  consequently,  the  spin  coherence  time  (and  length)  will  be  longer  in  our 
molecule.  An  order  of  magnitude  increase  in  lifetime  will  increase  the  coherence  length  by  a 
factor  ofVlO,  to  about  150  mn. 


Third,  the  work  function  of  MEFI-PPV  (5.1  eV)  is  nearly  well  matched  to  that  of  contact  LSMO 
(4.8  eV),  and  the  contact  resistance  in  our  device  is  expected  to  be  much  lower.  We  will  revisit 
these  calculations  when  data  on  our  structures  are  available.  To  study  the  trends,  we  use  two 
values  of  Ls  and  the  parameters  deduced  from  FET  measurements. 


From  FET  measurements,  the  1/f  noise  current  spectral  density,  Si,  was  found  to  be  proportional 

2  2 

to  i  Id  ,  where  i  is  the  total  current  and  d  is  the  channel  thickness.  With  Co  as  the  proportionality 
constant,  we  have 


S;  = 


C0i2 

d2/ 


(4) 


16  2 

From  the  measured  values  of  /,  d,  and  Si,  we  obtain  a  value  of  10'  cm'  for  Co  at/=l  Hz.  Then, 
from  the  well  known  and  simple  relationship,  the  1/f  noise  magnetic  spectral  density  Sh  and 
sensitivity  r)  in  N  x  N  array  of  sensor  can  be  obtained. 
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In  our  transport  modeling  work,  as  described  in  the  previous  section,  we  found  that  the  largest 
(dV/dB)  is  possible  only  when  the  field  is  constant  in  the  device.  The  constant  field  condition  is 
satisfied  when  the  doping  density  is  at  least  5x10  cm'  .  In  our  performance  studies  reported 
here,  we  assumed  a  doping  density  of  1  x  10 19  cm'3  and  used  a  linear  model  for  transport. 

For  a  chosen  value  of  1000  nm  for  d,  10  kV/cm  for  electric  field  strength  E,  we  calculate  the 
potential  V  across  the  organics  and  the  bias  magnetic  field  at  which  the  (dV/dB)  is  largest.  These 
values  are  substituted  in  Eq.  (6)  for  a  100  x  100  array  to  obtain  sensitivity  at  1  Hz  as  a  function 
of  carrier  mobility,  p,  and  plotted  for  two  values  of  spin  coherence  lengths  as  shown  in 
Figure  10. 

We  see  that  detectivity  (solid  lines)  improves  with  spin  coherence  length  and  the  best  value  of 

~10  fT/VHz  is  predicted  when  p  ~  5  x  10'7  cm2/Vs  in  this  sensor  array  operating  at  room 

2 

temperature.  Also  plotted  in  Figure  10  is  the  required  power  density  (dotted  lines)  in  mW/cm  . 
Best  sensitivity  with  very  low  power  consumption  (~5  mW/cm  )  is  possible  in  these  devices. 

So  far,  we  have  used  a  value  for  Ls  that  is  independent  of  mobility.  Since  the  diffusion  constant 
depends  on  mobility,  it  is  possible  that  Ls  is  proportional  to  Vp.  Under  this  assumption,  the 
calculations  are  repeated  and  results  are  displayed  in  Figure  1 1.  At  very  low  mobilities,  the  spin 
coherence  length  decreases  and,  consequently,  the  detectivity  performance  deteriorates.  The 
minimum  in  detectivity  is  ~  50  TT/a/Hz,  but  at  a  slightly  higher  mobility  (~  1 0  "5  cm2/Vs).  The 

required  higher  mobility  results  in  lower  resistance  and  higher  current.  Correspondingly,  the 

2 

required  power  density  also  increases  to  nearly  100  mW/cm 

The  effects  of  the  electric  field  and  doping  density  are  shown  in  Figure  12a.  With  the  smaller 
reduction  in  E  to  8  KV/cm  (blue  lines),  the  power  density  gets  smaller.  However,  the  sensitivity 
performance  gets  correspondingly  worse.  Similarly,  by  decreasing  the  dopant  density  from  1019 
cm"  to  10  cm"  ,  the  conductivity  decreases  and  hence  the  power  density  also  decreases.  Note 
that  the  sensitivity  perfonnance  gets  worse  by  a  factor  of  2.  Finally,  we  carried  out  studies  to 
evaluate  the  thickness  dependence.  The  results  are  shown  in  Figure  12b. 
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Figure  1 0:  Predicted  variation  of  detectivity  and  power  density  with 
mobility  and  spin  coherence  length. 


Figure  1 1 :  Predicted  variation  of  detectivity  and  power  density  with 
mobility  and  modified  spin  coherence  length. 


16 


Detectivity  r\  [pT/Hz 


Figure  12:  Predicted  variation  of  detectivity  and  power  density  with  electric  field,  dopant  density 
(a)  and  device  thickness  (b).  The  assumed  parameters  are  shown  in  the  legend. 


When  the  thickness  is  much  smaller  than  the  spin  diffusion  length,  we  do  not  take  full  advantage 
of  long  transit  times  required  for  enhanced  performance.  As  the  thickness  is  increased,  we  have 
more  time  to  manipulate  the  electron  spin  and  the  detectivity  improves.  Beyond  some  threshold 
value  of  thickness,  the  detectivity  saturates  at  about  60  LT/a/Hz.  Since  the  increase  in  length  also 
causes  the  resistance  to  increase,  the  power  density  increases  as  shown.  If  the  thickness  is  larger 
than  the  spin  diffusion  length,  both  detectivity  and  power  dissipation  will  be  worse. 

A  clarification  about  the  power  density  is  in  order.  Although  the  calculated  power  density  values 
appear  to  be  large  (~  100  mW/cm2),  this  heat  is  generated  only  across  the  organic  whose  cross- 
sectional  area  is  200  pm  x  0.3  pm.  So,  total  power  used  by  a  100  x  100  array  is  (100  x  10'  x  200 

8  4 

x  0.3  x  10'  x  10  )  W.  Since  this  array  is  expected  to  be  placed  on  2.5  cm  x  2.5  cm  substrates,  the 
fully  normalized  power  density  is  only  96  pW/cm“. 

1.9  Conclusions 

The  sensitivity  and  power  density  values  depend  intricately  on  the  device  thickness,  spin 
coherence  length,  mobility,  electric  field,  dopant  density,  and  bias  magnetic  field.  Although  it  is 
difficult  to  arrive  at  optimum  values  for  all  these  parameters  without  a  noise  model  appropriate 
to  our  polymer  and  device  design,  we  can  observe  the  following  requirements. 

•  The  surface  roughness  of  the  LSMO  layer  needs  to  be  improved.  The  preferred  target 
level  for  the  boulder  density  is  ~  104  per  cm2.  Since  we  plan  to  have  104  devices  on 
1  in.  x  1  in.  substrate,  we  would  on  the  average  have  one  boulder  per  LSMO  contact. 

•  With  improved  growth  quality  of  LSMO,  the  Tc  is  expected  to  slightly  increase  from  the 
current  value  of  360  K.  The  theoretical  maximum  is  370  K. 
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•  Currently  LSMO  is  grown  on  SrTi03  substrates,  which  are  commercially  available  only 
in  1  cm  x  1  cm.  However,  for  eventual  fabrication  of  a  100  x  100  array  on  a  1  in.  x  1  in. 
area,  LSMO  growth  on  large  area  substrates  must  be  explored. 

•  The  work  function  of  LSMO  is  known  to  vary  from  4.8  eV.  We  require  organic  materials 
whose  HOMO  level  (for  hole  devices)  or  LUMO  level  (for  electron  devices)  is  located 
near  5.1  eV.  This  will  reduce  the  contact  resistance  and  enable  efficient  charge  injection. 
The  HOMO  level  value  available  in  the  literature  for  MEH-PPV  is  5.1  eV  (i.e.,  0.3  eV 
barrier  for  holes).  With  appropriate  treatment  of  the  interface,  this  barrier  could  be 
overcome. 

•  The  exact  value  of  the  required  mobility  will  depend  on  the  noise  model  developed,  but  it 

7  5  2* 

appears  that  mobilities  in  the  range  5x10'  to  5  x  10"  cnr/Vs  will  be  required. 

•  Appropriate  approaches  to  reduce  the  contact  barrier  (by  about  0.5  eV)  and  improve 
ohmic  contact  must  be  explored.  For  example,  interface  dipoles  can  be  exploited. 

•  The  organic  must  be  truly  disordered  and  amorphous  in  form.  Appropriate  treatment  to 
either  remove  the  defects  or  bury  their  levels  deep  in  the  band  gap  must  be  explored. 
Charged  impurity  scattering  and  low-lying  defect  levels  contribute  to  noise. 

18  3 

•  The  n-doping  density  in  the  organic  must  be  at  least  5  x  10  cm’  ,  and  preferably 
10  cm'  .  The  electric  field  does  not  appear  to  be  an  issue,  because  dielectric  breakdown 
of  most  organics  occurs  at  ~  5  xlO6  V/cm,  and  our  device  must  operate  in  the  ~  104  V/cm 
range. 

•  Although  the  FIB  method  is  successful  in  making  required  nanotrenches,  the  use  of  Ga+ 
ion  negates  the  purpose  by  coating  the  LSMO  surface  with  oxide.  The  spin  and  charge 
injections  are  severely  affected. 

•  The  e-beam  lithography  with  Ar+  ion  milling  should  be  used  for  making  nano  trenches. 

•  Signal  analysis  techniques  for  an  efficient  on-chip  integration  of  100  x  100  array  signal 
must  be  developed. 

•  Since  the  response  function  takes  its  maximum  value  at  very  low  bias  magnetic  fields 
(~pG  to  mG),  this  array  of  devices  will  most  likely  be  operated  on  a  feedback  loop 
control  mode.  Explicitly,  an  external  constant  current  circuit  will  provide  the  required 
optimum  bias  field  in  the  organics.  When  the  signal  magnetic  field  is  detected,  the  current 
in  the  circuit  will  change  and  the  external  circuit  will  adjust  the  current  to  compensate 
this  change.  The  signal  is  deduced  from  the  correction  applied  by  the  external  circuit. 

2.  FIELD  EFFECT  TRANSISTOR 

Both  the  magnetic  sensor  and  the  magnetic-field-effect  transistor  have  a  common  LSMO- 
polymer-LSMO  structure.  Once  a  metal  and  a  semiconductor  are  brought  into  close  contact,  a 
Schottky  barrier,  as  well  as  a  depletion  layer,  will  be  formed  at  the  interface,  which  can 
significantly  modify  charge  transport  and  electric-field  distribution  inside  a  device.  Since  our 
device’s  operation  sensitively  depends  on  the  interplay  of  the  transit  and  precession  times,  it  is 
important  to  examine  the  robustness  of  the  performance  predicted  with  simplified  models. 
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We  have  developed  a  complete  device  model,  where  both  the  Schottky-barrier  effect  and  the 
electric-field-dcpcndcnt  carrier  mobility  are  taken  into  account.  We  have  considered  both  the 
thermionic-emission  and  the  tunneling  currents  and  solved  charge  and  spin  transport  equations 
self-consistently.  This  model  provides  a  realistic  and  accurate  description  of  spin-dependent 
transport  in  the  organic  spin  devices. 

2.1  Schottky  Barrier 

The  Schottky  barrier  height  in  a  metal-semiconductor  junction  is 

<l>Bn  =  <Kn-X 


where  e<f>m  is  the  metal  work  function,  and  ex  is  the  electron  affinity  measured  from  the  bottom 
of  the  conduction  band.  Due  to  the  image  force,  the  Schottky  barrier  is  lowered, 


$m  =  4>m 


where  E(0)  is  the  field  at  the  boundary  and  the  barrier  reduction  only  occurs  when  the  direction 
of  E(0)  is  pointing  from  the  metal  to  the  semiconductor.  The  current  across  a  barrier  consists  of  a 
contribution  from  thermionic  emission  and  electron  tunneling,  which  leads  to  a  spin-dependent 
interfacial  conductance.  The  up-spin  conductance  is 
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where  the  first  tenn  is  from  the  thermionic-emission  and  the  second  is  from  the  tunneling.  The 
Richardson  constant  A  is  given  by 

.  *  47rem*k2 


where  m  is  the  electron  effective  mass  of  the  polymer  and 

_e/z  (n7 
00  2  v  em* 

measures  the  effectiveness  of  tunneling,  which  depends  on  both  the  doping  concentration  Nd  and 
m*.  The  total  conductance  at  the  LSMO-polymer  interface  is  then 

G.=(l/Gj  +  1/Gt°j'  , 

where  G°  is  the  up-spin  conductance  in  the  absence  of  the  Schottky  barrier.  The  conductance  for 
different  spins  are  related  by  GJG'=  G’{'  /  Gf1  andGtyGt  =  /G® .  The  carrier  density  at  the 
boundary  between  the  metal  and  the  semiconductor  is 
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where  Nc  is  the  effective  density  of  state  in  the  semiconductor,  which  is  proportional  to  (in*)  . 
m*  in  a  polymer  is  usually  larger  than  in  an  inorganic  semiconductor,  such  as  Si. 

2.2  Electric-Field  Dependent  Mobility 

Carrier  mobilities  in  polymers  often  strongly  depend  on  the  electric  field  and  are  found  to  obey 
the  Poole-Frenkel  relation, 


where  v  is  the  mobility  and  Ec  measures  how  strong  the  field-dependence  is.  This  Poole-Frenkel 
behavior  originates  from  the  disorders  (traps)  in  the  polymer.  In  a  polymer  with  many  carriers,  as 
in  these  organic  spin  devices,  the  deep  traps  will  likely  be  filled  by  the  carriers,  and  the  field- 
dependence  is  not  as  strong  as  in  undoped  or  very  lightly  doped  polymers. 


2.3  Transport  Equations 

The  field-dependent  carrier  mobility  will  introduce  a  few  new  terms  in  the  charge-spin  transport 
equations.  We  solve  the  following  equations 

d2p0  dp0  f  e E  1  dv  dE2!  ep0  ( dE  E  dv  dE3 
dx‘  dx  ^kBT  v  dE  dx  J  kgT'ydx  v  dE  dx  J 

d2/  d/feE  ldvdE^  dln(p0)df  r  ?  /  n 

dx‘  dx  ^kBT  v  dE  dx  J  dx  dx  L2 

T= - (Po-Nd) 

dx  s 


where  /  =  p  /  ( epo/kBT ).  Here  both  the  spin  diffusion  length  L  =  -\[Dt as  well  as  the  effective 
magnetic  field  h  =  gjuBB / Dfi  depend  on  the  electric  field,  via  the  Einstein’s  relation  between  the 
mobility  and  the  diffusion  constant,  and  consequently  are  inhomogeneous  across  the  polymer. 
The  Schottky  effects  are  reflected  in  the  boundary  conditions,  as  described  earlier.  We  have 
developed  an  efficient  and  reliable  code  to  numerically  solve  these  equations  self-consistently. 


2.4  Results 

We  choose  parameters  pertinent  to  a  realistic  device  structure.  Nc  is  1021  cm'3,  which  is 
considerably  higher  than  in  Si  (3.8x10  cm'  ).  The  corresponding  carrier  effective  mass  of  the 
polymer  is  m*= 3.55  mo,  where  mo  is  the  bare  electron  mass.  The  Schottky  barrier  of  0.25  V  is 
assumed.  This  is  a  typical  energy  difference  between  work  functions  of  LSMO  and  MEH-PPV. 
A  high  barrier  would  not  allow  even  a  decent  charge  injection. 


20 


Figure  13  illustrates  the  electric-field  and  carrier-density  distributions  of  a  LSMO-polymer- 
LSMO  device.  With  a  value  of  Ec  set  1010  V/cm,  the  mobility  is  essentially  independent  of  the 
field.  We  see  that  a  depletion  layer  is  formed  at  both  interfaces  and  that  the  electric  field  in  the 
depletion  layers  significantly  deviates  from  the  value  in  the  middle  of  the  device.  Note  that  the 
two  depletion  layers  are  not  symmetric  because  in  the  sandwiched  structure  considered  here,  one 
Schottky  barrier  is  under  a  forward  bias  and  the  other  is  under  a  reverse  bias. 


Figure  13:  Spatial  distributions  of  electric  field  (upper  panel)  and  carrier  density  (lower 
panel)  for  an  LSMO-polymer-LSMO  device  structure.  The  Schottky  barrier 
height  is  0.25  V  and  the  doping  in  the  polymer  is  1018  cm"3. 


We  plot  in  Figure  14  the  device  sensitivity  dVdB  for  different  doping  concentrations  No.  For  a 
chosen  barrier  height  of  0.25  V,  we  see  that  the  dV/dB  value  obtained  in  jVd=1018  cm""1  case  is 
larger  than  that  obtained  in  the  case  of  either  No=  1 0 1 7  or  10 19  cm""1.  Since  the  average  electric 
field  as  well  as  other  parameters  is  fixed  in  these  calculations  for  these  doping  levels,  a  larger 
dV/dB  value  indicates  a  higher  spin  injection  efficiency  from  LSMO  into  the  polymer.  For  a 
slightly  different  barrier  height  of  0.3  V,  note  from  Figure  15  that  the  dV/dB  value  is  larger  for 
jVd=1017  cm"3  than  for  jVd=1018  or  1019  cm'3. 

We  plot  in  Figure  14  the  device  sensitivity  dVdB  for  different  doping  concentrations  No-  For  a 

i  o  -3 

chosen  barrier  height  of  0.25  V,  we  see  that  the  dV/dB  value  obtained  in  jVd=10  cm'  case  is 
larger  than  that  obtained  in  the  case  of  either  No=  1 0 1 7  or  10 19  cm""’.  Since  the  average  electric 
field  as  well  as  other  parameters  is  fixed  in  these  calculations  for  these  doping  levels,  a  larger 
dV/dB  value  indicates  a  higher  spin  injection  efficiency  from  LSMO  into  the  polymer.  For  a 
slightly  different  barrier  height  of  0.3  V,  note  from  Figure  15  that  the  dV/dB  value  is  larger  for 
jVd=1017  cm""1  than  for  jVd=1018  or  1019  cm""’. 
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Figure  14.  Device  sensitivity  dV/dB  (upper  panel)  and  device  voltage  V  (lower  panel)  as  a  function  of  an 
applied  magnetic  field  for  different  doping  concentrations.  Constant  current  and  electric  field 
(-5x104  V/cm,)  are  assumed.  The  barrier  height  is  0.25  V. 


Figure  15.  Device  sensitivity  dV/dB  (upper  panel)  and  device  voltage  V  (lower  panel)  as  a  function  of 
applied  magnetic  field  for  different  doping  concentrations.  The  barrier  height  is  0.3  V. 
Other  parameters  are  the  same  as  in  Fig.  14. 
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As  Nd  increases,  the  Fermi  level  in  the  polymer  moves  closer  to  the  conduction  band  edge.  The 
effective  barrier  height  increases  and  consequently  the  thermionic-emission  current  is  reduced 
while  the  tunneling  current  is  enhanced,  which  modifies  dV/dB.  Thus  the  above  results  indicate 
that  spin  injection  can  be  enhanced  by  engineering  Schottky  barrier  properties  such  as  the  barrier 
height  and  the  doping  in  the  semiconductor. 

Figure  16  compares  results  with  different  Ac  for  a  barrier  height  of  0.25  V.  The  values  dV/dB  for 
Ac  =  5x10"  and  10  cm'  are  similar,  but  are  quite  different  from  those  for  Ac  =  2x10  cm'  .  A 
smaller  Ac  means  a  smaller  carrier  effective  mass,  which  would  result  in  a  larger  tunneling 
current.  The  reduced  dV/dB  for  Ac  =  2x10  cm"  suggests  that  when  effective  electron  tunneling 
(smaller  effective  mass)  is  larger,  the  spin  injection  is  reduced. 

We  describe  the  effect  of  the  field-dependent  mobility  on  the  device  performance  in  Figure  17. 
We  see  that  this  effect  becomes  notable  only  when  Ec  is  105  V/cm,  and  even  then  dV/dB  is  not 
significantly  affected.  Ec=  105  V/cm  is  a  typical  value  for  undoped  polymers  and  denotes  a  very 
strong  electric-field  dependence  of  the  mobility.  It  appears  that  field  dependence  in  the  cases  of 
doped  polymers  will  not  play  a  significant  role.  The  weak  effect  of  field-dependent  mobility  on 
the  device  performance  is  possibly  because  the  electric  field  is  largely  homogeneous  inside  the 
device,  except  in  the  depletion  layers. 


Figure  16:  Device  sensitivity  dV/dB  (upper  panel)  and  device  voltage  V  (lower  panel)  as  a  function  of 

applied  magnetic  field  for  values  of  Nc.  The  barrier  height  is  0.25  V.  Other  parameters  are  the 
same  as  in  Figure  8. 
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Figure  17:  Device  sensitivity  dV/dB  (upper  panel)  and  device  voltage  V  (lower  panel)  as  a 
function  of  applied  magnetic  field  for  different  Ec,  which  measures  how 
sensitive  the  mobility  depends  on  the  field.  The  barrier  height  is  0.25  V.  Other 
parameters  are  the  same  as  in  Figure  14. 


In  conclusion,  our  results,  which  are  based  on  the  realistic  and  accurate  modeling  that  includes 
Schottky  barrier  effects  as  well  as  field  dependence  of  the  mobility,  have  demonstrated  the 
robustness  of  the  performance  predicted  from  a  simplified  model.  This  model  also  suggests  new 
degrees  of  freedom,  such  as  the  barrier  height,  carrier  effective  mass,  and  doping  concentrations, 
to  be  engineered  to  achieve  high  performance  of  organic  spin  devices. 

3.  LIGHT  PROPAGATION  IN  GaAs 

SRI  had  previously  developed  codes  to  study  high-intensity  light  propagation  through 
semiconductor  compounds  and  alloys.  The  codes  are  now  extended  to  the  short- wavelength 
region  and  applied  to  GaAs.  The  required  propagation  parameters,  including  intrinsic  carrier 
density,  one-photon  absorption  coefficient  (a),  two-photon  absorption  coefficient  (P),  free  carrier 
absorption  coefficient  (a),  and  refractive  index  changes  arising  from  carrier  density  and 
temperature,  are  calculated  from  full  band  structures.  The  excellent  comparison  between  the 
calculated  and  the  measured  a  (as  seen  in  Figure  18)  indicates  that  both  band  structures  and 
dipole  matrix  element  calculations  are  accurate.  The  calculated  values  of  p  are  small,  of  the  order 
of  10  cm/GW,  and  decrease  with  the  wavelength  (Figure  19). 
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Figure  18:  One-photon  absorption  in  GaAs 


Figure  19:  Two-photon  absorption  in  GaAs 


Due  to  the  absence  of  free  carrier  absorption  by  holes  in  the  wavelength  region  of  1  pm  to  2  pm, 
the  calculated  a  is  extremely  small  (of  the  order  of  10'  cm  ).  However,  the  cross  section 
increases  with  wavelength  and  temperature  (shown  in  steps  of  50  K  in  Figure  20).  The  calculated 
values  of  a,  P,  a  are  used  in  a  Kramer-Kronig  relationship  to  obtain  the  change  in  refractive 
index.  We  find  that  the  index  changes  linearly  with  the  carrier  density  (Figure  21)  with  a  slope  of 
-2  x  10'21.  The  calculated  refractive  index  increases  with  T  at  a  rate  of  about  10‘4  per  K. 


Wavelength  [gm] 

Figure  20:  Free  carrier  absorption  in  GaAs 


Figure  21 :  Index  change  with  carrier  density  in  GaAs 
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These  calculated  values  are  then  used  in  our  three-dimensional  light  propagation  code  to  obtain 
the  output  intensity  as  a  function  of  input  intensity.  This  code  considers  both  absorption  and 
refraction  in  the  material  on  numerical  mesh  and  solves  the  time-dependent  Maxwell’s  equations 
as  a  function  of  distance.  The  calculated  limiting  of  the  light  with  a  wavelength  of  1.06  pm  is 
shown  in  Figure  22  (solid  line);  the  measured  values  [Boggess  et  al.,  IEEE  J.  Quan.  Elec.  21,  488 
(1985)]  are  shown  as  filled  circles.  We  see  that  the  trend  is  predicted  accurately  and  the 
calculated  values  are  within  experimental  accuracy.  If  the  [1  value  of  22  cm/GW  is  used  instead 
of  18  cm/GW,  we  see  (dashed  line)  an  excellent  numerical  agreement  as  well. 


Input  Intensity  [GW/cm2] 


Figure  22:  Optical  limiting  in  GaAs. 


The  codes  can  now  be  used  to  study  1.5  pm  light  propagation  through  organic  molecules  and 
polymers  with  the  values  of  [:>  and  a  measured  under  DARPA’s  MORPH  program.  In  the  next 
period,  we  will  evaluate  the  nonlinear  absorption  parameters  required  for  the  optical  limiting 
goals  of  the  MORPH  program.  We  have  further  applied  our  codes  to  study  light  propagation 
through  GalnAs  and  GaAsP  alloys  and  concluded  that  GaAsP  alloys  are  better  than  either  GaAs 
or  GalnAs  alloys  for  limiting  1.5  pm.  The  results  are  summarized  for  possible  publication  in 
Applied  Physics  Letters.  The  draft  manuscript  is  enclosed  here  as  Appendix  B. 
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We  present  a  comprehensive  model  to  describe  spin  transport  across  an  organic  material  sand¬ 
wiched  between  magnetic  contacts.  We  solve  a  set  of  differential  equations  for  transport  of  charge 
and  spin,  together  with  the  Poisson  equation.  Both  the  electric-field-induced  spin  and  charge  drift 
and  the  magnetic-field-induced  spin  precession  are  included.  We  also  consider  the  effects  of  a 
Schottky  barrier  between  the  metal  contact  and  organic  materials  and  of  electric-field  dependent 
carrier  mobility.  The  magnet-polymer-magnet  structures  were  proposed  to  achieve  ultrasensitive 
magnetometers  and  low-power  magnetic-field-effect  transistors  because  a  weak  magnetic  field  can 
significantly  modify  spin  transport  through  spin  precession.  The  numerical  results  based  on  the 
comprehensive  modeling  suggest  that  the  organic  device  structures  that  exploit  low  carrier  mobility 
and  long  spin  lifetimes  for  magnetic-sensing  and  magnetic-field-effect  transistors  are  robust. 

PACS  numbers: 


I.  INTRODUCTION 

In  spin-based  electronics  (spintronics),  information  is 
carried  by  the  electron  spin.  Because  of  the  nature 
of  magnetic  interactions  between  electron  spins,  spin- 
tronic  devices  have  the  potential  advantages  of  non¬ 
volatility,  increased  data  processing  speed,  and  decreased 
electric  power  consumption.  Inorganic  semiconductor 
spintronic  devices  have  attracted  considerable  attention1 
since  the  discovery  of  long  spin  lifetimes  in  semiconduc¬ 
tor  structures2.  Injection,  manipulation,  and  detection 
of  spin  polarization  in  semiconductors  have  been  exten¬ 
sively  investigated  both  experimentally  and  theoretically. 

Compared  with  inorganic  materials,  organics  have 
much  longer  spin  lifetimes  because  of  the  van¬ 
ishing  spin-orbit  couplings,  suggesting  that  organic 
materials  have  significant  potential  for  novel  spin¬ 
tronic  devices.  Recently,  strong  magnetoresistances 
and  large  spin  injection  have  been  observed  in 
Lao.7Sro.3MnC>3  (LSMO)/sexithienyl  (Tg)/LSMO  and 
LSMO/8-hydroxyquinolate  aluminum  (Alq3)/Co  struc¬ 
tures  even  at  room  temperature3,4.  T6  and  Alq3  are  two 
widely  used  materials  in  organic  electronics  like  light- 
emitting  diodes.  Theoretical  studies  of  spin-dependent 
transport  in  magnet-polymer-magnet  structures  have 
just  begun.  A  systematic  theoretical  study  of  spin  in¬ 
jection,  spin  manipulation,  and  spin  detection  in  organic 
structures  is  required  to  understand  the  experiments  and 
to  design  new  organic  spintronic  devices. 

Organic  electronic  devices,  including  light-emitting 
diodes  and  field-effect  transistors,  have  been  the  sub¬ 
jects  of  intense  research  in  the  last  decade  because  they 
have  processing  and  performance  advantages  for  low-cost 
and/or  large-area  applications5.  On  the  other  hand,  the 
low  carrier  mobility  (1CU8  —  10  cm2/Vs)  in  organics  lim¬ 
its  their  application  for  high-speed  devices  such  as  com¬ 
puter  processors.  The  arena  of  organic  electronic  devices 
is  in  the  applications  that  exploit  unique  properties  in  or¬ 
ganics,  such  as  large-area  processing,  mechanical  flexibil¬ 
ity,  tunable  light  emission,  chemical  sensing  interactions, 


and  biocompatibility.  Similarly,  organic  spintronic  de¬ 
vices  have  some  unique  applications  that  are  inaccessible 
to  inorganic  spintronic  devices. 

Recently  a  spin  drift-diffusion-precession  equation  was 
derived  to  explain  experiments  and  design  organic  spin 
devices6.  It  was  found  that  the  low  carrier  mobility  and 
the  long  spin  lifetime  in  organic  materials  provide  an 
ideal  condition  to  exploit  the  interplay  of  electric-field- 
induced  spin  drift  and  magnetic-field-induced  spin  pre¬ 
cession  to  make  ultrasensitive  magnetic  sensors  and  low- 
energy  magnetic-field-effect  transistors.  Very  recently, 
spin  precession  and  the  resulting  oscillatory  resistance 
were  successfully  observed  in  the  spin-drift  regime7  in 
Si8.  In  Ref.  6,  to  capture  major  physics,  a  few  approx¬ 
imations  that  assume  a  uniform  electric  field  across  the 
device  structure  and  a  carrier  mobility  independent  of 
electric  field  were  adopted.  These  approximations  allow 
one  to  obtain  separate  linear  differential  equations  for 
spin  and  charge  transport,  which  are  relatively  easy  to 
solve.  In  general,  however,  the  electric  field  is  inhomo¬ 
geneous,  particularly  near  the  contacts  in  organic  device 
structures.  In  addition,  carrier  mobility  in  organic  ma¬ 
terials  may  strongly  depend  on  electric  field  and  carrier 
density.  Moreover,  there  exists  a  Schottky  barrier  be¬ 
tween  the  metal  contact  and  the  organic  material.  Since 
the  devices  operation  exploits  the  interplay  of  the  tran¬ 
sit  and  precession  times,  it  is  important  to  examine  the 
robustness  of  the  performance  predicted  with  simplified 
models.  The  comprehensive  model  also  provides  a  de¬ 
signing  tool  for  organic  spintronic  devices. 

In  this  work  we  go  beyond  the  approximations  and 
provide  a  complete  description  of  spin-dependent  trans¬ 
port  in  magnet-polymer-magnet  structures.  To  make  the 
theory  general  enough  to  cover  most  structures,  we  in¬ 
clude  both  electric-field-induced  spin  drift  and  magnetic- 
field-induced  precession  and  the  magnetization  of  the  two 
magnets  can  be  noncollinear.  The  comprehensive  model 
also  provides  a  designing  tool  for  organic  spintronic  de¬ 
vices.  We  assume  spin-dependent  transport  across  or¬ 
ganic  structures  is  diffusive  rather  than  tunneling.  The 
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FIG.  1:  Spatial  distributions  of  electric  field  (upper  panel) 
and  carrier  density  (lower  panel)  for  an  LSMO-polymer- 
LSMO  device  structure.  The  Schottky  barrier  height  is  0.25 
V  and  the  doping  in  the  polymer  is  1018  cm-3. 
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FIG.  2:  Device  sensitivity  dV/dB  (upper  panel)  and  device 
voltage  V  (lower  panel)  as  a  function  of  an  applied  magnetic 
field  for  different  doping  concentrations.  Constant  current 
and  electric  field  (~  5  x  104  V/cm,)  are  assumed.  The  barrier 
height  is  0.25  V. 


paper  is  organized  as  follows.  In  Sec.  II,  a  complete  set 
of  drift-diffusion-precession  equations  are  presented.  In 
Sec.  Ill,  we  analyze  the  effects  of  the  Schottky  barrier 
and  the  field-dependent  mobility.  We  present  numerical 
results  in  Sec.  IV,  and  Sec.  V  summarizes  our  results. 


II.  SPIN  CHARGE 

DRIFT-DIFFUSION-PRECESSION  EQUATIONS 

When  a  voltage  is  applied  to  a  magnet-polymer- 
magnet  structure,  a  spin-polarized  current  is  injected 
into  the  polymer  from  the  magnets,  giving  rise  to  spin 
accumulation  in  the  polymer.  To  consider  spin  preces¬ 
sion  and  noncollinear  configurations,  where  spin  accu¬ 
mulation  can  be  along  any  direction,  we  use  a  2x2  den¬ 
sity  matrix  in  spin  space  to  describe  the  carrier  distri¬ 
bution,  pp  =  ppl  +  &  ■  pF .  Here  pp  1  is  the  equilib¬ 
rium  carrier  distribution  of  the  nonmagnetic  polymer, 
and  <r  =  (<fx, <ry, az)  are  Pauli  matrices.  In  the  previ¬ 
ous  work,  to  emphasize  the  spin-dependent  part  in  car¬ 
rier  transport,  as  a  simplification,  we  assume  that  the 
charge  distribution  inside  the  polymer  is  homogeneous, 
and  V pp  =  0  and  V  ■  E  =  0,  and  in  the  steady  state  we 
obtain 


where  b  =  gpBB/HD  and  L  =  y/Drs-  This  equation 
provides  a  simple  description  of  spin  drift  and  spin  pre¬ 
cession  in  polymers. 

In  general,  however,  the  local  electric  field  and  the 
charge  distribution  are  interdependent,  an  accurate  de¬ 
scription  requires  self-consistently  solving  Poisson’s  equa¬ 
tion  together  with  transport  equations.  To  this  end,  we 
start  from  the  general  form  of  a  spin-polarized  current  in 
the  organic  material,  which  consists  of  two  contributions, 


drift  and  diffusion, 

jp  =  ppevE  -  eDX7pp,  (2) 

where  v  is  the  carrier  mobility  and  D  the  diffusion  con¬ 
stant  in  the  organic  material.  Here  we  neglect  the  pos¬ 
sible  magnetic-field  effect  on  the  orbital  motion  (Hall  ef¬ 
fect),  which  is  reasonable  in  polymers  with  low  carrier 
mobilities.  In  a  nondegenerate  system,  v  and  D  are  con¬ 
nected  via  Einstein’s  relation  v/eD  =  l/kBT.  The  conti¬ 
nuity  equation  for  each  component  of  the  density  matrix 
in  the  presence  of  a  magnetic  field,  B,  reads 

dpp  pp-ppi  1  .p  i 

dt  ts  e  h 


99B 


(a-B)  ,  (3) 


The  above  transport  equations  need  to  be  solved  to¬ 
gether  with  the  Poisson  equation, 

X7  .  E  = -^(p0  ~  Nd)  (4) 

where  NB  is  the  ionized  donor  density.  The  inclusion  of 
the  Poisson  equation  in  the  spin  transport  theory  makes 
an  exact  solution  impossible,  and  we  have  to  resort  to 
numerical  integration  of  these  coupled  differential  equa¬ 
tions.  Suppose  the  electric  field  is  along  the  axaxis  and 
spin  and  charge  distributions  are  uniform  along  y  —  z 
plane,  the  spin-charge  drift-diffusion  equations  in  the 

rs  -P 

steady  state  (gp-  =  0)  can  be  written  as 

d2po  dpo  /  eE  1  dv  dE\  epo  /  dE  E  dv  dE\ 
dx2  dx  \kBT  +  v  dE  dx  /  kBT  V  dx  v  dE  dx  ) 

(5) 


d2p 

dx2 


/  eE  1  dv  dE\ 

1  |  eP  I 

rdE 

E  dv  dE\ 

V  kBT  1  v  dE  dx  ) 

+  kBT ' 

y  dx 

v  dE  dx  ) 
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Dh 


B  x  p 


(6) 
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FIG.  3:  Device  sensitivity  dV/dB  (upper  panel)  and  device 
voltage  V  (lower  panel)  as  a  function  of  applied  magnetic  field 
for  different  doping  concentrations.  The  barrier  height  is  0.3 
V.  Other  parameters  are  the  same  as  in  Fig.  2. 


FIG.  4:  Device  sensitivity  dV/dB  (upper  panel)  and  device 
voltage  V  (lower  panel)  as  a  function  of  applied  magnetic 
field  for  values  of  Nc ■  The  barrier  height  is  0.25  V.  Other 
parameters  are  the  same  as  in  Fig.  2. 


Here  the  carrier  mobility  is  assumed  to  be  a  function  of 
electric  field,  v  =  v(E).  Introducing  chemical  potentials 
f  for  spins,  p  =  Eq.  (6)  can  be  rewritten  as 

d2f  df  /  eE  1  dv  dE\  d\np0  df  f  = 

dx 2  dx\kBT  v  dE  dx )  dx  dx  L 2 

(7) 

where  both  h  and  L  are  electric-field  dependent  through 
v  and  D,  h  =  gpsB /hD(E)  and  L  =  y J D{E)ts ■ 
Equations  (5)  and  (7)  represent  the  interdependent 
spin-charge  transport  equations,  which  need  to  be  solved 
self-consistently  with  the  Poisson’s  equation  and  proper 
boundary  conditions. 


III.  ELECTRIC-FIELD  DEPENDENT 
MOBILITY 

Carrier  mobilities  in  polymers  often  strongly  depend 
on  the  electric  field  and  are  found  to  obey  a  universal 
Poole-Frenkel  relation, 


v  =  Uq  exp 


(8) 


where  v  is  the  mobility  and  Ec  measures  how  strong  the 
field-dependence  is.  This  Poole-Frenkel  behavior  origi¬ 
nates  from  the  disorders  (traps)  in  the  polymer.  In  a 
polymer  with  many  carriers,  as  in  these  organic  spin  de¬ 
vices,  the  deep  traps  will  likely  be  filled  by  the  carriers, 
and  the  field-dependence  is  not  so  strong  as  in  undoped 
or  very  lightly  doped  polymers.  Using  the  Poole-Frenkel 
relation,  we  have 


E  dv  _  1  l~E 
v  dE  2  V  Ec 


IV.  SCHOTTKY  BARRIER  AND  BOUNDARY 
CONDITIONS 


The  two  magnets  in  a  magnet-polymer-magnet  can  be 
regarded  as  magnetic  reservoirs  in  local  equilibrium  at 
chemical  potentials  Pen'  which  is  diagonal  in  spin  space 

Pen  =  Pen)-  Here  C(TZ)  denotes  the  left  (right)  mag¬ 
net.  The  direction  of  the  magnetization  in  each  magnet 
is  described  by  the  unit  vector  me,n-  In  the  absence  of 
Schottky  barrier  effects  the  current  leaving  the  left  con¬ 
tact  at  x  =  0  and  the  current  entering  the  right  contact 
at  x  =  d  are9,10 


JC(0) 


GWc\fic  -  AP(°)]4  +  Gl0ulc[pe  -  AP(0)]«£ 

Gllulpp(0)ulc  -  GlJulcpp(0)ul  (10) 


jC(d )  =  Glu]n[pp(d)  -  Pn]un  +  G^[pp(d)  -  p^n 
+  G^u^id^  +  G^u^id)^.  (11) 

Here  pp  is  the  polymer  electrochemical  potential  in  the 
spin  space.  Operators  ^  \  [l  +  (— )<r  •  mc\  and 

|  [l  +  (— )<r  ■  mn]  project  spins  to  the  magne¬ 
tization  directions  of  the  magnets.  The  above  equations 
can  be  regarded  as  a  generalized  Ohm’s  law  in  the  spin 
space.  Gj  (Go)  is  the  electron  conductance  in  the  mag¬ 
net  with  spin  parallel  (antiparallel)  to  the  magnetization 
direction.  Gj^  =  ReGj^  +  ilmGy'  is  the  mixing  conduc¬ 
tance,  which  measures  the  transport  capability  of  spins 
oriented  perpendicular  to  the  magnetization  direction. 

In  a  magnet-polymer-magnet  structure,  generally 
speaking,  a  Schottky  barrier,  as  well  as  a  depletion  layer, 
will  be  formed  at  the  interface  between  the  magnet  and 
polymer,  which  can  significantly  modify  charge  transport 
and  electric- field  distribution  inside  a  device.  The  Schot¬ 
tky  barrier  height  in  a  metal-semiconductor  junction  is 

4*Bn  =  0m  X  (12) 


which  will  be  used  when  solving  Eqs.  (5)  and  (7). 


where  ecf>m  is  metal  work  function,  and  e\  is  the  elec¬ 
tron  affinity  measured  from  the  bottom  of  the  conduc¬ 
tion  band.  Due  to  the  image  force,  the  Schottky  barrier 
is  lowered, 


where  -E(O)  is  the  electric  field  at  the  boundary.  The 
barrier  reduction  only  occurs  when  the  direction  of  E( 0) 
is  pointing  from  the  metal  to  the  semiconductor.  The 
current  across  a  barrier  consists  of  a  contribution  from 
thermionic  emission  and  electron  tunneling,  which  leads 
to  a  spin-dependent  interfacial  conductance.  The  up-spin 
conductance  is 


gt  =  A*Ti 


Bn 

~k^f 


a&Bn  \ 

EqO  )  -  ksT 


(14) 


where  the  first  term  is  from  the  thermionic-emission  and 
the  second  is  from  the  tunneling. 
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FIG.  5:  Device  sensitivity  dV/dB  (upper  panel)  and  device 
voltage  V  (lower  panel)  as  a  function  of  applied  magnetic  field 
for  different  l?c,  which  measures  how  sensitive  the  mobility 
depends  on  the  field.  The  barrier  height  is  0.25  V.  Other 
parameters  are  the  same  as  in  Fig.  2. 
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The  Richardson  constant  A*  is  given  by 

47 Tem*k2 


A*  = 


h 3 


(15) 


(16) 


where  m*  is  the  electron  effective  mass  of  the  polymer 
and  measures  the  effectiveness  of  tunneling,  which  de¬ 
pends  on  both  the  doping  concentration  Nd  and  m* . 
The  total  conductance  at  the  magnet-polymer  interface 
is  then 


GT  =  (1/Gj +  1/GI)-1.  (17) 

If  the  difference  between  Gj  and  Gq  is  due  mainly  to 
the  energy  gap  between  up-spin  and  down-spin  electrons 
in  the  magnet,  as  in  most  half  metals,  the  relations 
for  conductance  Go  with  difference  spins  would  remain 
valid  for  G  with  different  spins,  G^/G^-  =  GJ/Gq  and 

gti/gt  =  g^/gJ. 

Inside  the  polymer,  the  currents  of  spin  and  charge  can 
be  expressed  as 


jo 
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ep0vE  +  vksT- 
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e2povE  dp0  df 

— ~  veETJ  ~  euPo~r- 
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(18) 

(19) 


Substituting  the  four  spin-dependent  components  of  Go 
in  Eqs.  (8)  and  (9)  by  those  of  G,  we  have  obtained  a 
set  of  new  boundary  conditions  for  the  magnet-polymer- 
magnet  structure  with  a  Schottky  barrier. 

The  condition  of  a  continuous  current  is  not  sufficient 
to  determine  the  charge  distribution  in  the  depletion  re¬ 
gion  of  a  Schottky  barrier.  One  must  specify  the  charge 


density  of  the  polymer  at  the  at  the  boundary  between 
the  magnet  and  the  polymer, 


k^T 


eE(  0) 
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where  E(  0)  is  the  electric  field  at  x  =  0,  which  reduces 
the  barrier  height  through  the  image  charge,  pc  is  the 
effective  density  of  state  in  the  semiconductor,  which 
is  proportional  to  (m*)3/2.  m*  in  a  polymer  is  usually 
larger  than  in  an  inorganic  semiconductor,  such  as  Si. 


V.  RESULTS 

We  have  developed  an  efficient  and  accurate  numeri¬ 
cal  approach  to  solve  the  coupled  transport  equations  of 
spin  and  charge  together  with  the  Poisson  equation.  The 
magnetic  contact  is  assumed  to  be  LSMO,  a  half-metal, 
with  a  ration  of  G^/G-*-  =  107.  We  choose  parameters 
pertinent  to  a  realistic  device  structure.  pc  =  1021  cm-3, 
which  is  considerably  higher  than  in  Si  (3.8  x  1019  cm-3). 
The  corresponding  carrier  effective  mass  of  the  polymer 
is  m*  =  3.55too,  where  mo  is  the  bare  electron  mass.  The 
Schottky  barrier  of  0.25  V  is  assumed.  This  is  a  typical 
energy  difference  between  work  functions  of  LSMO  and 
MEH-PPV.  From  our  calculation,  a  high  barrier  height 
larger  than  0.5  V,  would  not  allow  even  a  decent  charge 
injection  from  LSMO  into  a  polymer. 

Figure  1  illustrates  the  electric- field  and  carrier-density 
distributions  of  a  LSMO-polymer-LSMO  device.  With  a 
value  of  Ec  set  1010  V/cm,  the  mobility  is  essentially  in¬ 
dependent  of  the  field.  We  see  that  a  depletion  layer  is 
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formed  at  both  interfaces  and  the  electric  field  in  the  de¬ 
pletion  layers  significantly  deviates  from  the  value  in  the 
middle  of  the  device.  Note  that  the  two  depletion  layers 
are  not  symmetric  because  in  the  sandwiched  structure 
considered  here  one  Schottky  barrier  is  under  a  forward 
bias  and  the  other  is  under  a  reverse  bias. 

We  plot  in  Figure  2  the  device  sensitivity  dV/dB  for 
different  doping  concentrations  Nd-  For  a  chosen  barrier 
height  of  0.25  V,  we  see  that  the  dV/dB  value  obtained  in 
Nd  =  1018  cm-3  case  is  larger  than  that  obtained  in  the 
case  of  either  Nd  =  1017  or  1019  cm-3.  Since  the  average 
electric  field  as  well  as  other  parameters  is  fixed  in  these 
calculations  for  these  doping  levels,  a  larger  dV/dB  value 
indicates  a  higher  spin  injection  efficiency  from  LSMO 
into  the  polymer.  For  a  slightly  different  barrier  height 
of  0.3  V,  note  from  Figure  3  that  the  dV/dB  value  is 
larger  for  Nd  =  1017  cm-3  than  for  Nd  =  1018  or  1019 
cm-3.  As  Nd  increases,  the  Fermi  level  in  the  polymer 
moves  closer  to  the  conduction  band  edge.  The  effective 
barrier  height  increases  and  consequently  the  thermionic- 
emission  current  is  reduced  while  the  tunneling  current  is 
enhanced,  which  modifies  dV/dB.  Thus  the  above  results 
indicate  that  spin  injection  can  be  enhanced  by  engineer¬ 
ing  Schottky  barrier  properties  such  as  the  barrier  height 
and  the  doping  in  the  semiconductor. 

Figure  4  compares  results  with  different  Nq  for  a  bar¬ 
rier  height  of  0.25  V.  The  values  dV/dB  for  Nc  =  5  x  1021 
and  1021  cm-3  are  similar,  but  are  quite  different  from 
those  for  Nc  =  2  x  102°  cm-3.  A  smaller  Nc  means 
a  smaller  carrier  effective  mass,  which  would  result  in 
a  larger  tunneling  current.  The  reduced  dV/dB  for 
Nc  =  2  x  102°  cm-3  suggests  that  when  effective  elec¬ 
tron  tunneling  (smaller  effective  mass)  is  larger,  the  spin 
injection  is  reduced.  We  describe  the  effect  of  the  field- 
dependent  mobility  on  the  device  performance  in  Figure 
5.  We  see  that  this  effect  becomes  notable  only  when  Ec 
is  105  V/cm,  and  even  then  dV/dB  is  not  significantly 
affected.  Ec  =  105  V/cm  is  a  typical  value  for  undoped 
polymers  and  denotes  a  very  strong  electric-field  depen¬ 


dence  of  the  mobility.  It  appears  that  field  dependence 
in  the  cases  of  doped  polymers  will  not  play  a  significant 
role.  The  weak  effect  of  field-dependent  mobility  on  the 
device  performance  is  possibly  because  the  electric  field 
is  largely  homogeneous  inside  the  device  except  in  the 
depletion  layers. 

In  conclusion,  our  results,  which  are  based  on  the  real¬ 
istic  and  accurate  modeling  that  includes  Schottky  bar¬ 
rier  effects  as  well  as  field  dependence  of  the  mobility, 
have  demonstrated  the  robustness  of  the  performance 
predicted  from  a  simplified  model.  This  model  also  sug¬ 
gests  new  degrees  of  freedom,  such  as  the  barrier  height, 
carrier  effective  mass,  and  doping  concentrations,  to  be 
engineered  to  achieve  high  performance  of  organic  spin 
devices. 

VI.  CONCLUSIONS 

We  have  presented  a  comprehensive  theory  to  describe 
spin  transport  in  magnet-polymer-magnet  structures. 
This  theory  considers  both  the  electric-field-induced  spin 
drift  and  magnetic-field-induced  spin  precession.  Our  re¬ 
sults,  which  are  based  on  the  realistic  and  accurate  mod¬ 
eling  that  includes  Schottky  barrier  effects  as  well  as  field 
dependence  of  the  mobility,  have  demonstrated  the  ro¬ 
bustness  of  the  performance  predicted  from  a  simplified 
model.  This  model  also  suggests  new  degrees  of  free¬ 
dom,  such  as  the  barrier  height,  carrier  effective  mass, 
and  doping  concentrations,  to  be  engineered  to  achieve 
high  performance  of  organic  spin  devices. 
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We  have  developed  a  code  to  study  nonlinear  propagation  of  high  intensity  light  through  semi¬ 
conductor  compounds  and  alloys.  Nonlinear  absorption  and  refraction  are  included.  The  input 
parameters-  absorption  and  refraction  coefficients,  temperature  and  carrier  densities-  are  inter¬ 
dependent  and  change  in  the  sample  as  the  light  propagates.  Using  full  bandstructures  of  GaAs,  we 
evaluated  all  required  parameters  for  1.06  /trn  wavelength  light.  The  transmission  ratio  of  the  output 
intensity  to  the  input  intensity,  calculated  without  any  adjustable  parameter,  agree  very  well  with 
measured  values.  We  further  calculated  the  parameters  for  InP,  Ga64lii36As,  and  I11P60AS40  alloys 
to  study  the  propagation  of  1.55  /im  wavelength  light  and  predict  largest  attenuation  of  transmitted 
intensity  in  Ga64lii36As  alloy. 

PACS  numbers:  78.67.ch,42.65.k 


Studies  of  nonlinear  absorption  (NLA)  and  nonlinear 
refraction  (NLR)  in  semiconductors  continues  to  be  of 
scientific  and  technological  interest.1,2  NLA  is  used  in  a 
number  of  applications,  including  frequency  conversion, 
optical  switches  and  optical  limiting.1  Consequently,  the 
computer  codes  to  model  light  propagation  through  the 
medium  have  been  developed  specifically  for  the  chosen 
application. 3-12  The  modeling  of  output  intensity  and 
beam  profile  variation  with  input  intensity  in  semicon¬ 
ducting  materials  is  being  carried  out  for  scientific  as 
well  as  the  selection  of  better  materials  for  optical  limit¬ 
ing  applications.4,5  In  these  models,  the  NLA  and  NLR 
are  characterized  by  two-photon  absorption  (TPA)  and 
free  carrier  absorption  (FCA)  coefficients,  and  nonlin¬ 
ear  refraction  coefficients.  However,  in  reality  the  val¬ 
ues  of  these  coefficients  change  as  the  light  propagates 
through  the  material,  owing  to  the  creation  of  large  num¬ 
ber  of  photo-carriers  and  generation  of  heat.  For  exam¬ 
ple,  when  more  carriers  are  generated,  more  band  edge 
states  are  filled  causing  a  drastic  reduction  in  the  TPA 
coefficients.13,14  Similarly,  the  changes  in  both  carrier 
density  and  lattice  temperature  affect  the  magnitude  and 
sign  of  the  change  in  the  refractive  index. 15-18  The  com¬ 
putational  method  and  code  to  study  the  beam  propaga¬ 
tion  should  include  the  transient  nature  of  all  nonlinear 
parameters. 

In  this  Letter,  we  (a)  outline  the  method  used  in  our 
studies  of  beam  propagation  through  nonlinear  medium 
in  three-dimensions,  (b)  use  all  inter-dependent  param¬ 
eters  calculated  from  full  bandstructures  of  GaAs,  (c) 
demonstrate  excellent  agreement  between  the  calculated 
transmission  and  measured  values,  and  (d)  apply  this 
method  to  study  the  propagation  of  1.55  /im  light  beam 
through  GaAs,  InP,  Ga64ln36As,  and  InP60As4o  alloys. 
While  the  limiting  is  larger  in  smaller  gap  materials,  we 


predict  higher  performance  in  Ga64ln3gAs  alloy  than 
that  in  in  InPgoAs4o  alloy  of  nearly  equal  band  gap. 

The  details  of  the  calculational  method  for  propagation 
will  be  published  elsewhere.19  Here  we  give  a  brief  out¬ 
line  for  continuity.  Since  the  study  of  three  dimensional 
propagation  of  light  requires  both  magnitude  and  phase 
of  the  electric  field,  we  start  from  the  form  of  Maxwell’s 
equation,  instead  of  the  rate  equation,  and  include  non¬ 
linear  polarization  term.20 


V2E(r,f)  =  ~n 


d2E(Y,t)  dE(r,t)  d2PNL{r,t) 


82t  c  dt 
PNL(r,t)  =  3e0x(3)(w,w,  -aj)\E(r,t)\2E(r,t) 


d2t 


(1) 


We  assume  p=po,  and  e=n2e o,  where  /to  and  eo  as  the 
permeability  and  permittivity  of  the  free  space  respec¬ 
tively  and  n  is  the  refractive  index  of  the  medium.  In 
the  presence  of  FCA,  the  conductivity  ac  changes  to 
<rc o  +  Acrc  and  n  changes  to  no  +  An,  where  the  sub¬ 
script  0  denote  the  corresponding  value  in  the  unper¬ 
turbed  medium.  We  then  substitute  E(r,t)  =  A(r)  e_ifcz 
e~lut  in  Eq.  1  and  use  the  relationships  between  the 
one  photon  absorption  coefficient  a  and  a,  the  nonlinear 
refractive  index  7  and  the  real  part  of  the  third-order 
susceptibility  x^3\  and  the  TPA  coefficient  (3  and  the 
imaginary  part  of  x ^  to  get 
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The  subscript  T  denotes  the  transverse-to-propagation 
direction  and  the  corresponding  component  in  the  prop- 
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FIG.  1:  Radial  distribution  of  the  beam  just  inside  the  entry 
surface  (dotted  line)  and  exit  surface  with  (thick  solid  line) 
and  without  (thin  solid  line)  refraction  terms.  Input  pulse  is 
normalized  to  1  at  the  the  beam  center 

agation  direction  is  assumed  to  be  zero  (paraxial  ap¬ 
proximation).  af  is  the  FCA  coefficient,  and  ko  is  the 
wavevector  in  the  unperturbed  medium.  The  first  term 
determines  the  magnitude,  whereas  the  second  and  third 
terms  determine  the  phase  of  the  electric  field  of  the  prop¬ 
agating  wave.  Eq.  2  is  solved  on  a  three-dimensional 
numerical  mesh  by  rewriting  it  in  an  implicit  form  and 
care  is  taken  to  achieve  numerical  stability  by  using 
staggered-leapfrog  differencing  scheme.  With  measured 
incident  pulse  profile  as  the  initial  boundary  condition, 
the  beam  profile  at  the  exit  surface  is  calculated  as  a 
function  of  time  and  radial  distance  by  solving  Eq.2  self- 
consistently  with  the  rate  equations13,14  governing  the 
photo-generated  carrier  density  (A N)  and  temperature 
(T)in  the  sample.  The  propagation  of  the  beam  profile 
from  the  exit  surface  to  the  aperture  is  calculated  in  Han- 
kel  function  approach21  and  integrated  as  appropriate  for 
closed  or  open  aperture  measurements. 

An  important  ingredient  in  our  beam  propagation 
modeling  is  the  use  of  (3,  a/ ,  and  An,  calculated 
from  full  bandstructures,  along  with  their  complicated 
dependence14  on  AN  and  T.  Since  AN  and  T  depends  on 
r  and  t,  these  parameters  are  continuously  evaluated  and 
used  in  the  equations  as  the  beam  propagates  through 
the  sample.  Our  previous  calculations  in  the  midwave 
infrared  (MWIR)  regime  indicates  strong  variation  of  the 
parameters  (3  and  FCA  cross  section  a  (which  is  the  ra¬ 
tio  of  FCA  coefficient  to  carrier  density)  with  T  and  AN. 
The  calculations  carried  out  here  in  the  1.06  /jin  to  1.55 
/an  regime  also  indicate  strong  variation  with  T,  but  not 
with  AN  for  densities  up  to  1018  cm-3. 

The  nonlinear  parameters  (/3,  tr,  and  An)  in  GaAs  are 
calculated  as  described  previously  from  the  underlying 
tight-binding  bandstructures.22  At  MWIR  wavelengths, 


FIG.  2:  The  transmission  calculated  without  any  adjustable 
parameter  (solid  line)  is  compared  to  the  measured  values 
(data).  The  value  of  3  adjusted  to  yield  good  fit  (dashed 
line)  is  22  cm/GW 


the  value  of  (3  and  cr  in  semiconductors  are  of  the  order  of 
1  cm/MW  and  10~16  cm2.  However,  the  corresponding 
values  at  1.06  /an  wavelength  in  intrinsic  GaAs  are  cal¬ 
culated  to  be  1.8  xlO-2  cm/MW  (or  18  cm/GW)  and  7.1 
x  ICC20  cm2  at  300K.  Since  the  material  is  required  to 
have  larger  band  gap  for  nonlinear  absorption  at  shorter 
wavelengths,  the  [3  is  very  small.  The  FCA  cross  section 
is  small  because  of  the  holes,  which  play  a  dominant  role 
in  FCA  at  MWIR  wavelengths,  can  not  participate  in  the 
absorption  of  short  wavelength  photons.  In  other  words, 
the  electrons  in  light-hole  valence  band  do  not  find  a  final 
state  at  the  same  wavevector  k  in  heavy  hole  band  after 
the  direct  absorption  of  1.2  eV  photons.  Both  electrons 
in  the  conduction  band  and  holes  in  the  valence  bands 
contribute  to  FCA  only  through  an  indirect  process,  with 
the  assistance  from  phonons,  and  hence  the  cross  section 
is  extremely  small. 

The  calculated  nonlinear  parameters  at  1.06  /im  are 
used  in  the  beam  propagation  model  to  study  the  trans¬ 
mission  of  laser  pulse  with  half- width  at  1/e  of  the  max¬ 
imum  (HWl/eM)  of  27  ps  and  60  /iin.  These  beam  pa¬ 
rameters,  the  sample  thickness  (of  0.175  cm)  and  inci¬ 
dent  intensity  values  are  taken  from  Ref.23.  The  sample 
is  not  anti-reflection  (AR)  coated.  A  measured24  value  of 
-4.1xl0-16  cm2/W  for  7  is  used.  The  incident  beam  just 
inside  the  GaAs  entry  surface  is  shown  in  Fig.  1  by  dot¬ 
ted  line.  Also  shown  are  the  calculated  beam  profile  at 
the  exit  surface  with  (thick  solid  line)  and  without  (thin 
solid  line)  the  refraction  terms.  Since  the  absorption  is 
proportional  to  intensity,  more  of  the  light  is  absorbed 
at  the  center  of  the  beam.  Hence  we  see  that  even  when 
the  refraction  is  not  considered,  the  transmitted  beam  is 
non-gaussian.  In  the  presence  of  refraction,  the  beam 
spreads  additionally  and  away  from  the  center,  resulting 
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TABLE  I:  Calculated  nonlinear  parameters  in  intrinsic  mate¬ 
rial  at  T=300K  and  A=1 .55/un 


Parameters  GaAs  InP  InP(joAs4o  Ga64lii36As 


(3  [cm/GW] 

7.3 

12.3 

32.9 

36.2 

a  [10-20  cm2] 

21.6 

14.7 

23.0 

22.1 

dn/dN  [10~21cm3] 

-1.05 

-2.3 

-15.9 

-16.2 

dn/dT  [10-4KC1] 

1.3 

0.54 

2.1 

2.5 

in  highly  non-gaussian  beam  profile  as  shown. 

The  transmitted  output  energy  is  simply  the  time  and 
radial  integration  of  beam  profile  at  the  exit  surface.  The 
transmission  coefficient,  which  is  the  ratio  of  output  en¬ 
ergy  to  input  energy,  calculated  with  full  bandstructure 
parameters  is  plotted  (thick  solid  line)  in  Fig.  2  as  a 
function  of  input  intensities.  Both  trend  and  magnitude 
compare  well  with  the  measured  values  (filled  circles).23 
However,  when  the  (3  value  of  18  cm/GW  at  300K  is 
slightly  increased  to  22  cm/GW,  we  can  obtain  a  near 
perfect  fit  (dashed  line)  with  experimental  values.  Hence 
we  conclude  that  both  calculated  and  experimental  val¬ 
ues  agree  very  well  within  the  numerical  and  measure¬ 
ment  accuracies. 

The  calculations  are  then  carried  out  to  predict  materi¬ 
als’  limiting  performance  at  A=1.55  /iin.  For  the  compar¬ 
ative  study,  we  considered  GaAs  and  InP  compounds  and 
Ga64ln36As  and  InPgoAs4o  alloys.  GaAs  and  InP  have 
near  equal  bandgap  of  1.42  eV  and  1.35  eV  respectively 
at  300K.  The  alloy  concentrations  in  Ga64ln3gAs  and 
InP60As40  are  chosen  so  that  they  have  equal  band  gap 
(of  0.92  eV)  at  300K.  The  nonlinear  parameters  /3,  cr, 
dn/dN ,  and  dn/dT  obtained  from  underlying  bandstruc- 
tures  at  300K  are  given  in  Table  I.  We  note  that  the 
P  in  GaAs  is  reduced  when  the  wavelength  is  increased 
from  1.06  to  1.55  ;im,  and  that  it  increases  when  the 
bandgap  is  decreased.  Although  the  bandgap  in  InP  is 
nearly  equal  to  that  in  GaAs,  the  calculated  /3  differs 
by  as  much  as  50%.  In  spite  of  near-identical  bandgap 
value  in  the  alloys,  the  /?  differ  by  about  10%.  The  FCA 
cross  section  is  known25  to  vary  as  A3  and  we  see  that 
the  calculated  a  in  GaAs  follows  that  variation.  We  fur¬ 
ther  find  that  a  is  nearly  independent  of  the  bandgap  of 
materials  considered.  As  before,  their  variation  with  T 
and  AN  are  also  calculated  and  used  in  our  code.  For 
simplicity,  the  incident  beam  profile  is  assumed  to  be 
Gaussian  with  temporal  and  spatial  HWl/eW  of  27  ps 
and  60  /mi  respectively,  the  sample  is  not  AR  coated 
and  has  a  thickness  of  1.75  mm  (all  parameters,  except 
wavelength,  are  identical  to  those  used  to  obtain  Fig.  1 
and  2).  The  transmission  coefficients  obtained  as  func¬ 
tions  of  incident  intensity  in  GaAs,  InP,  Ga64ln36As  and 


InP6oAs4o  are  shown  in  Fig.  3.  Because  of  the  difference 
in  the  refractive  index(  of  3.31  in  GaAs  and  3.1  in  InP) 
and  the  resulting  Fresnel  reflection,  the  transmission  at 
low  intensities  in  GaAs  and  InP  are  different.  At  higher 
intensities,  InP  is  expected  to  limit  better  than  GaAs.  In 
general,  we  see  that  the  limiting  at  a  given  wavelength 
is  enhanced  with  the  reduction  in  the  bandgap.  Con¬ 
siderably  enhanced  limiting  performance  is  predicted  in 
InP6oAs4o  and  Ga64ln36As  alloys.  Owing  to  slightly 
larger  /3,  the  Ga64ln36As  alloy  is  predicted  to  offer  the 
best  optical  limiting  performance  at  this  wavelength.  We 
call  for  experimental  verification  of  this  prediction. 


Input  Intensity  [GW/cm2] 

FIG.  3:  Calculated  transmission  of  1.55  /jm  beam  through 
GaAs,  InP,  Ga64lii36As  and  InP6oAs4o  as  a  function  of  in¬ 
cident  intensity. 

In  summary,  we  have  developed  a  method  to  study 
the  beam  propagation  in  three  dimension  by  including 
nonlinear  absorption  and  refraction.  When  the  inter¬ 
dependent  nonlinear  parameters  obtained  from  full  band- 
structures  are  used  in  the  code,  an  excellent  agreement 
with  measured  transmission  of  1.06  /mi  beam  through 
GaAs.  The  further  calculations  carried  out  for  1.55 
/mi  beam  predict  a  superior  limiting  performance  in 
InP60As40  and  Ga64ln36As  alloys,  with  GalnAs  offering 
a  slightly  better  performance. 
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